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A recent work 1 demonstrated, for an ideal spinless p+ip superconductor, that Majorana end-states 
can be realized outside the strict one-dimensional limit, so long as: 1) the sample width does not 
greatly exceed the superconducting coherence length and 2) an odd number of transverse sub-bands 
are occupied. Here we extend this analysis to the case of an effective p+ip superconductor engineered 
from Rashba spin-orbit coupled surface with induced magnetization and superconductivity, and find 
a number of new features. Specifically, we find that finite size quantization allows Majorana end- 
states even when the chemical potential is outside of the induced Zeeman gap where the bulk material 
would not be topological. This is relevant to proposals utilizing semiconducting quantum wires, 
however, we also find that the bulk energy gap is substantially reduced if the induced magnetization 
is too large. We next consider a slightly different geometry, and show that Majorana end-states can 
be created at the ends of ferromagnetic domains. Finally, we consider the case of meandering edges 
and find, surprisingly, that the existence of well-defined transverse sub-bands is not necessary for 
the formation of robust Majorana end-states. 

PACS numbers: 71.10.Pm, 74.20.Rp, 74.78.-w, 03.67.Lx 



I. INTRODUCTION 

Majorana fermion bound states are expected to ex- 
hibit non-Abelian exchange statistics 2 ' 3 , and have been 
proposed as a basis for topological quantum computers 
which would be protected from decoherence 4,5 . Conse- 
quently, there is a growing interest in producing Majo- 
rana fermions in the laboratory. Superconductors with 
p+ip pairing symmetry have long been expected to posses 
zero-energy Majorana bound states in vortex cores 3 . 
Such p + ip superconductors are thought to naturally oc- 
cur in triplet paired fermionic superfluids (such as 3 He A 
or Sr2Ru04) 6 ' 7 , and in the Pfafflan quantum Hall state at 
v = 5/2 8 . However, despite extensive experimental work 
on such systems, direct evidence of Majorana fermions 
remains elusive. 

Recently, the possibility of engineering effective p + 
ip superconductors in more conventional materials has 
arisen 9 12 . A particularly promising class of such pro- 
posals involves using Rashba-type spin-orbit coupling in 
combination with conventional s-wave superconductivity 
to produce an effective p x + ip y 2D superconductor 9 11 . 
Magnetization would then be introduced to remove one 
of the two components, leaving an effective p x + ip y su- 
perconductor. 

The practical difficulties of creating and manipulat- 
ing vortices has led to a renewed interest in the original 
Kitaev idea 5 , where Majorana particles are realized as 
localized states at the ends of a one-dimensional p x + ip y 
superconducting wire. In addition to being potentially 
simpler to implement than vortex based proposals, cre- 
ating Majorana fermions as end-states in quantum wires 
would allow one to build scalable networks of gates to 
braid, fuse, and measure many Majorana fermions 11 . A 
recent work 1 demonstrated, for an ideal (spinless) p + ip 
superconductor, that Majorana end-states can be real- 
ized outside the strict one-dimensional limit, so long as: 



1) the sample width does not greatly exceed the super- 
conducting coherence length £o = 7tuf/A, and 2) an odd 
number of transverse sub-bands are occupied. Further- 
more, since Majorana end- states emerge only for an odd 
number of occupied sub-bands, as chemical potential is 
swept the system undergoes a sequence of alternating 
topological phase transitions between phases with and 
without Majorana end-states. 

In this paper we extend this analysis to the case of 
spinfull fermions with Rashba spin-orbit coupling and 
with s-wave pairing and ferromagnetic splitting induced 
by proximity effect. Because this system is effectively a 
single species p + ip superconductor for a certain range of 
parameters, we expect, and indeed find, that the results 
of Ref. 1 still apply. However, because of the presence of 
multiple energy scales (Fermi-energy £p, spin-orbit cou- 
pling Ar, magnetization V z , and superconductivity A), 
several new features and possibilities emerge. 

We begin by analyzing a long narrow strip geometry 
with hard- wall boundary conditions. This geometry was 
previously considered in Ref. 13, which identifies Majo- 
rana end-states by computing a topological invariant for 
a small number of occupied sub-bands. Here, we extend 
these results to an arbitrary number of subbands. We 
find that, in the parameter regime where the system is 
effectively a single component p+ip superconductor, Ma- 
jorana end-states exist when an odd number of transverse 
sub-bands are occupied. Interestingly for V z > A > A^, 
a regime that is relevant to semiconductor materials, fi- 
nite width quantization allows the Majorana end-states 
to persist even for /jl > V z where a two-dimensional sam- 
ple would be topologically trivial. This observation al- 
lows one to operate at substantially larger carrier den- 
sity than previously expected 11 , placing less stringent re- 
quirements on sample purity. However, our results also 
show that the energy gap is substantially reduced if the 
Zeeman energy gets too large. 
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We also show that spatially non-homogenous magne- 
tization profiles can be used to produce Majorana end- 
states at the ends of long rectangular Ferromagnetic do- 
mains. We further demonstrate that structure of Majo- 
rana end-states remains largely unchanged for smoothly 
varying magnetization profiles. Such geometries are ad- 
vantageous for microfabricated structures, because the 
induced magnetization profile will be smooth even if the 
etched edges of the ferromagnetic material are rough, 
thus diminishing the impact of edge disorder. 

In addition, a similar setup could be used to real- 
ize end-states in an island of topological insulator (TI) 
with induced superconductivity (SC), surrounded by fer- 
romagnetic (FM) insulator (see Fig. 1). One can take 
advantage of the sensitivity of the existence of Majo- 
rana end states to the chemical potential by adopting 
the geometry shown in Figure 1. In this geometry, the 
Majorana modes could be moved around by selectively 
applying gate voltages to locally tune the number of oc- 
cupied sub-bands, thus obviating the need to create and 
manipulate vortices. 

Finally, we consider random, meandering edge geome- 
tries in order to address the question of whether or not it 
is necessary to have well-defined transverse sub-bands in 
order to produce Majorana end-states. Suprisingly, we 
find that the existence of Majorana fermions and the al- 
ternating structure of topological phase transitions with 
chemical potential persists, even when edge variations 
are large enough that there is no well-defined concept of 
transverse sub-bands. 

The paper is organized as follows: we begin with a 
short review of the proposed route to engineering an ef- 
fective single species p+ip superconductor from materials 
with Rashba spin-orbit coupling, and of the results of Ref. 
1 for quasi-one-dimensional spinless p + ip superconduc- 
tors. We then introduce the tight-binding model which 
forms the basis of our analysis, and describe a Green's 
function based method used to treat large system-sizes. 
We then describe the resulting analysis of this tight- 
binding model for hard- wall boundary conditions, ferro- 
magnetic domains, and random non-rectangular samples. 



II. OVERVIEW: TOPOLOGICAL 
SUPERCONDUCTIVITY FROM RASHBA 
COUPLED SURFACES 
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FIG. 1. (Color online) Proposed setup for electrically manip- 
ulating Majorana end-states in topological insulator (TI) ma- 
terials. A strip of TI with induced superconductivity (labeled 
TI+SC) is embedded in a ferromagnetic insulator (labeled 
FM) . Top gates (shown as overlayed rectangles) are used to lo- 
cally control the number of occupied sub-bands. Blue shaded 
gates indicate an even number of sub-bands, demarking a non- 
topological region, whereas un-shaded gates indicate an odd 
number of sub-bands, demarking a topological region. Ma- 
jorana bound states (shown as yellow blobs) emerge at the 
boundary between topological and non-topological regions. 
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where 



creates two helical bands with energies s± 

aji\k\ and spin- wavefunct ions = ( ^ 

0k = tan -1 (k x /k y ), which wind counter-clockwise and 
clockwise respectively. One can introduce s-wave (spin- 
singlet) superconductivity by proximity effect: 



kt C -kJ. 



h.c. 



(2) 



Re-expressing H/\ in terms of the helical Rashba surface 
bands, one finds that that the induced pairing has p±ip 
for respectively 9 ' 10 . 

While the Rashba splitting has effectively generated p- 
wave superconductivity, fermion states on the spin-orbit 
coupled surface still occur in degenerate time-reversed 
pairs. Consequently, the system is still topologically triv- 
ial and will not exhibit Majorana states. To break this 
time-reversal doubling and obtain a topologically non- 
trivial single-species p + ip superconductor, one can in- 
troduce Zeeman splitting term 



v z 



kf C kt 



(3) 



In this section we briefly review the proposed route 
to engineering an effective p + ip superconductor in 2D 
surfaces with Rashba spin-orbit coupling 9,10 . The Hamil- 
tonian for a 2D surface with Rashba coupling is: 



H R . 



ashba 



k 



[6c + QLRZ • (cr X k)] a0 4 a C k/3 (1) 



7 2 

where = ^ — fi is the spin-independent band struc- 
ture, ji is the chemical potential, an is the Rashba cou- 
pling strength, a are the spin-1/2 Pauli matrices, and 
a,/? =tl are spin indices. The Rashba coupling cxr 



through, for example, proximity to a Ferromagnetic in- 
sulator. Such a term modifies the bare (A = 0) energy 
bands to 



e± = 6 ± \ 



a 2 R k 2 



(4) 



V z also tends to cant the helical bands out of the xy- 
plane, giving them some component along the z-axis. 
Re-expressing in the eigenbasis of both Rashba and 
Zeeman couplings, one finds that, in addition to p db ip 
pairing A p (/c)k ± ~ (ck,±C-k,±) between fermions both in 
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band e±, this canting introduces an s-wave pairing com- 
ponent A s (fc) ~ (ck,+C-k,-) between fermions c + and c_ 
in bands £ + and E- respectively where: 

/A s (k)\ = 1 ( V z \ ( ) 

V A p( /c )y 2y/V?+a 2 R k 2 V" a W v J 

and k = (k y ± zfe^) /fc. 

As discussed in 11 , one has a topological superconduc- 
tor with potential Majorana bound states so long as V z > 
A, and so long as \i lies within the Zeeman gap < V z ). 
The latter restriction is potentially problematic for real- 
izing the above outlined scheme in semiconductor het- 
erostructures. These structures exhibit small Rashba 
splittings on the order of aRkp = 2ma 2 R ~ 10 _4 eV 
(where kp is the Fermi momentum for \i = 0). Fur- 
thermore, it is also desireable to have Rasbha splitting 
comparable or larger than Zeeman splitting such that 
the induced superconductivity has a substantial p-wave 
component (see Eq.5). For very small Rasbha splitting, 
the conditions < V z and aRkp > V z together require 
low carrier density, making such structures susceptible to 
disorder. 



III. MAJORANA END-MODES IN RASHBA 
COUPLED STRUCTURES 

A. Spinless p + ip Case 

Before treating the case of spinfull fermions with 
Rashba spin-orbit coupling, which will be the focus of 
this paper, we briefly review some pertinant results from 
Ref. 1 for Majorana end-states in spinless p+ip super- 
conductors. Ref. 1 considers a spinless p+ip supercon- 
ducting sample of length L x and width L y . In the 2D 
limit where L Xl L y ^> £o (where £o = ttvf/A is the su- 
perconducting coherence length), the sample has a bulk 
superconducting gap but exhibits a gapless chiral edge 
mode whose energy is quantized by the finite sample 
perimeter. The edge mode wave- function is localized on 
the system boundary with characteristic length-scale £o- 
As L y is decreased below £o> the tails of the edge mode 
wave- functions begin to strongly overlap generating a gap 
along the length of the sample that scales as ~ e~ Ly ^°. 

In Ref. 1, it was shown that for an odd number of oc- 
cupied transverse sub-bands, this Abuik ~ e~ Ly ^° bulk 
gap stabilizes zero-energy Majorana states isolated on 
opposite ends of the sample. Because these end-states 
occur when an odd number of transverse sub-bands are 
occupied, as chemical potential \i is changed, the sys- 
tem undergoes alternating sequence of topological phase 
transitions between phases with and without Majorana 
end-states. 

When Majorana end states are spatially well sepa- 
rated, the bulk gap Abuik exponentially suppresses the 
probability of an unpaired electron tunneling between 
any two Majorana states. This exponential suppression 



protects the Fermion parity information that is stored 
non- locally between any two Majoranas throughout any 
braiding process in which the particles remain well sep- 
arated. However, in order to measure the mutual occu- 
pation of two Majorana end-states, it is usefull to bring 
the Majoranas close together, fusing them into a single 
fermion state which is either occupied or unoccupied. For 
example, one can measure the occupation of two Majo- 
ranas by coupling them accross a Josephson junction and 
measuring the sign of the resulting Josephson current 5 . 
When fusing Majoranas for measurement, the important 
energy scale is not the gap to extended bulk-excitations 
but rather the so-called mini-gap to localized Fermion 
end-states. 

When Majoranas are realized as bound states in 
Abrikosov vortices, this mini-gap scales as Amg- Vortex ~ 
A 2 /sf <C A. For simple measurement schemes, this 
small mini-gap requires working at very low temper- 
ature. More sophisticated measurement schemes that 
do not require bringing Majoranas close to each other 
are possible 14 , however these schemes are comparatively 
more complex, presenting additional challenges for exper- 
iment. In comparison, the mini-gap for Majorana end- 
states in quasi-one-dimensional wires scales as 

A„o ~ *jj (6) 

where M is the number of occupied sub-bands. The slow 
square-root dependence on the number of occupied sub- 
bands indicates that the mini-gap can be a sizeable frac- 
tion of the bulk gap even for many occupied sub-bands. 
Figure 2 demonstrates this ^= for exact diagonalization 
of a quasi-one-dimensional spinless p + ip superconduc- 
tor. Panel A and B of Figure 2 show the scaling of the 
mini-gap with A and M respectively. 

The parametric scaling of the mini-gap can be under- 
stood simply, by considering each sub-band as contribut- 
ing a Majorana end-state, which overlap spatially, and 
are coupled to one another by the p-wave pairing gap 
~ A. Specifically, the A/\0\f excitation scaling is consis- 
tent with that of a random AfxAf antisymmetric matrix, 
whose entries are normally distributed with width A. 

In what follows, we demonstrate that a similar picture 
holds for effective p + ip superconductors generated from 
spinfull fermions with Rashba spin-orbit coupling, Zee- 
man splitting, and induced s-wave pairing. However, the 
story is complicated because of the presence of multiple 
energy scales (Fermi energy, Rashba and Zeeman split- 
tings, and superconducting gap) which gives rise to new 
features, and allows for a number of generalizations. 

B. Tight-Binding Model and Green's Function 
Method 

To analyze the structure of Majorana zero-modes in 
p+ip superconductors engineered from Rasbha coupled 
structures, we study numerically a discrete square-lattice 



FIG. 2. Demonstration that the mini-gap Amg protecting 
the Majorana end-statesof a spinless p + ip superconductor 
scales as A/VJ7 rather than A 2 /sf as is the case for vortex 
core states. The simulations were perfomed for a L x xL y strip 
of the tight-binding model considered in Ref. 1 with hopping 
t and p-wave BCS pairing A. Panel A shows the linear scaling 
of Amg with A for L x = 150, L y = 10, t = 10 and chemical 
potential fi = — 2t (corresponding to J\f = 5 filled sub-bands). 
Panel B shows the scaling of the mini-gap (where J\T 

is the number of occupied subbands), for L x — 250, L y — 20, 
t = 40, and A = 1. In Panel B, the dashed line shows the 
best power-law fit to the simulation data scales as iV -0-501 . 



tight-binding model version of the continuum Hamilto- 
nian: 

Htb = H t + Hso + i^FM + Ha 

H t = ~ l ( c R+d,c c R,« + ft.c.) - ^ na c na 

R,d,a 



H FM = V z C k* ( a ^*(3 C */3 + h ' C ' 



H so = ~ ia R c n+d,aZ • x d ) c R,/5 + h.c. 

R,d,a:,/3 

ff S C = X) Ac Rt c IU + ft - c - ( 7 ) 
R 

where R labels lattice sites, d G {e Xl e y } is a unit vector 
connecting nearest neighboring sites, (a, /3) G {tl} are 
spin-indices, and {&j}j= x ,y,z are spin-1/2 Pauli matrices. 
These terms represent the kinetic hopping energy (H t ), 
induced ferromagnetic Zeeman splitting (#fm), Rashba 
spin-orbit coupling (iJso) 5 an d induced spin-singlet pair- 
ing (Ha) respectively. 

For sufficiently small system size, the full 2D tight- 
bonding model (7) can be analyzed directly by exact di- 
agonalization. For larger systems, where exact diagonal- 
ization becomes computationally intractable, we employ 
a recursive Green's function method 17 depicted schemat- 
ically in Fig. 3. Starting with the Green's function for 
the right (left) end of a strip of length n, G^) L ^(y, y') one 
can construct the Green's function for the right (left) end 
of a strip of length n + 1 by: 



w + ir} - H n - VtG^V 
1 

uj + ir)- H n -VG { L ] Vi 



FIG. 3. (Color online) Schematic depiction of recursive 
Green's function method (eq. 8,9). (a) Green's function for 
the end of an ±jx~ layer sample G^ Lx ^> is calculated recursively 
from G (n_1) G (n) by adding one transverse layer at a time, 
(b) Green's function for any layer can then be recovered from 
solving (9) for the desired location, using the intermediate 
results for q( l *~ x ) anc [ from the top procedure. 



where H n is the Hamiltonian for the strip at location 
x = n, V is the matrix containing all elements of (7) that 
connect x = ntox = n + l, and 77 > gives the poles of 
G a small imaginary component. By repeated recursion 
of (8), one can find the Green's function for the end of 
an arbitrarily long strip. The Green's function for any 
fixed value of x can then be calculated from: 



G {Lx) (x,x;y 1 ,y 2 ;uj) 



cu + iri-H x - VtG^ -1) V - VG 



(L x — x) 



(9) 



(8) 



Once calculated, the Green's function can yield the den- 
sity of states: p(uj) = — ^Tr^mG(o;), which has a 5- 
function peak at each energy level. Furthermore, for 
uj — )> £ n , where e n is an eigen-energy of (7), the Green's 
function becomes a projector onto the corresponding 
eigenstate * n (R): lim^-^ (irj)G(R, R; uj) |^ n (R)| 2 . 



IV. HARD- WALL CONFINEMENT 

We first consider the the tight-binding Hamiltonian (7) 
on a long strip of length L x and width L y with hard- 
wall boundaries. These boundary conditions correspond 
to strong confinement and are relevant, for example, to 
self-assembled semiconductor nano- wires or etched semi- 
conductor or metallic microstrips. 

The bulk excitation gap, Abuik, for an infinitely long 
strip of width L y provides a convenient way of charac- 
terizing topological phase transitions. The locations of 
topological phase transitions, corresponding to closings 
of the bulk-gap, can be simply obtained by solving the 
effectively ID problem of diagonalizing (7) at zero mo- 
mentum along x (k x = 0). Using the Green's function 
method outlined above, we verify that closings of the bulk 
gap (Abuik = 0) in an infinitely long strip signify topo- 
logical phase transitions in a finite length strip. These 
transitions are between states with zero-energy Majo- 
rana end- modes, and topologically trivial states without 
zero-energy end-modes. Fig. 4 shows the topological 
phase-diagram as a function of chemical potential ji and 
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FIG. 4. (Color online) Phase diagram as a function of chem- 
ical potential fi and Zeeman splitting V z for strip of width 
L y — 10, with hard- wall boundaries and t = 10, A = 1, 
and q>r = 2. White regions correspond to the topologically 
trivial phase with no Majorana end-states. Colored regions 
indicate topologically non-trivial phases with Majorana end- 
states, color indicates size of bulk excitation gap Abuik which 
sets the length scale £m ~ i>f /Abuik over which the end- 
states decay along the length of the wire. Note the reduction 
of Abuik as V z becomes large. Dashed lines show the location 
of the non-superconducting (A = 0) transverse sub-bands, 
which are closely related to the location of topological phase- 
transitions. 



Zeeman splitting V z for a strip of width L y — 10 with 
parameters t = 10, A = 1, and cxr = 2. Fig. 4 is com- 
patible with the results of Ref. 13, which computes the 
phase diagram restricted to the lowest two sub-bands. 

The effect of transverse confinement in the y-direction 
is to quantize y- moment a k y to discrete values: k n with 
n G {1,2,...}, each with a corresponding a ID transverse 
sub-band. In Fig. 4, the transverse sub-band bottoms for 
the infinite strip with A = are shown as dashed lines. 
These transverse sub-bands are doubly spin-degenerate 
for V z = 0, and are split linearly as V z is increased. With- 
out Rashba splitting, this would result in a fan of crossing 
levels as V z is increased, however the presence of Rashba 
splitting leads to avoided crossings between the different 
sub-bands. 

Similarly to the case of multi-band of spinless p+ip su- 
perconducting quantum wires 1 , we find that zero-energy 
Majorana end-states exist only when an odd-number of 
these transverse sub-bands are occupied. However, the 
picture here is complicated by the presence of additional 
energy scales OLRkp and V z which give rise to two spin- 
split surface bands with energies s± = ± \/V z + a 2 R k 2 . 
As described above, pairing between two electrons of the 
same helicity has p x =p ip y symmetry and amplitude A p , 
whereas pairing between electrons of different helicity has 
s-wave symmetry with amplitude A s (see eq. 5). When 
V z is small enough that £+(fc n ) — S-(k n ) < A, both p + ip 
and p — ip pairing occurs and the system is topologically 
trivial for all //. In contrast, for sufficiently large V z such 
that e+(k n ) — E-{k n ) > A only one of the two p±ip com- 
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FIG. 5. Spectral function for the end layer of a 200 x 10 
strip for the tight-binding model (7) with hard-wall bound- 
aries as determined by Green's function method (8) with 
r\l A = 0.001. Simulation parameters used were: t = 10, 
A = 1, olr = 2, and V z = 6. From left to right and top to 
bottom, the chemical potentials are ji/t = -4.37, -4.09, -3.72, 
and -3.265 corresponding to Af = 1,2, 3 and 5 transverse sub- 
bands occupied respectively. There is a zero-energy pole only 
for Af odd. For general Af there are Af end-states split by the 
effective p-wave gap ~ A p . 



ponents remains. In this regime, so long as the quantized 
levels of e±(k y = k n ) are offset, there are intervals of fi 
for which an odd number, Af of occupied sub-bands. In 
such intervals there are pairs of p + ip and p — ip bands, 
and one unpaired p + ip band. One can picture each sub- 
band as contributing a single Majorana end-state; these 
are intercoupled and mix, to form /2)\ full fermions 
at non-zero energy and A/"mod2 Majorana zero- modes. 

By these considerations alone, it would seem that 
larger V z is always favorable for creating Majorana end- 
states. However, the gap protecting the Majorana end- 
states from bulk excitations depends on the p-wave com- 
ponent A p of the induced pairing gap, which for V z ^> 

2 

ma 2 R scales as: Abuik — A p — >> ^^A <C A (see Eq. 
5). Consequently, there is a trade off between increas- 
ing V z to stabilize Majorana end-states, and in avoiding 
V z ^> aR to protect the bulk pairing gap. 

This picture of each sub-band contributing a Majo- 
rana end-state is supported by Fig. 5a.-d. which show 
the spectral function for a L x = 200, L y = 10 strip, pro- 
jected onto the last layer. Fig. 5a.-d. were simulated us- 
ing the same parameters as for Fig. 4 with fixed V z = 6, 
and fi adjusted such that Af = 1, 2, 3, and 5 transverse 
sub-bands are occupied respectively. In each subfigure, 
one finds two types of excitations. The first are bulk ex- 
citations above Abuik ~ 0.25 A. These reside along the 
length of the strip and have small projection ~ 1/L X onto 
the end-layer. The second are end-states. The end state 
wave- functions are exponentially localized to the sam- 
ple ends, with localization length £m °° vf /Abuik set by 
the bulk excitation gap. Consequently, these states have 



6 



large projection onto the end-layer that is independent of 
L x for L x ^> £m- As shown in Fig. 5 for J\f occupied sub- 
bands, there are J\T end-states, and we have verified that 
they are roughly equally spaced in energy ~ A p /\^JV. 
For J\f odd, one of these end states sits at exactly zero 
energy, and is a Major ana state protected from excita- 
tions by energy gap: ~ A p . For Af even the end-states 
sit at non-zero energy and are not Major ana fermions. 

To summarize, we find that, in hard- wall confined 
strips, Majorana end states exist as long as the £ + — S- > 
A and so long as there are an odd number of occu- 
pied transverse sub-bands. Surprisingly, this alternating 
even-odd behavior persists in confined systems even for 
\fjb\ > V z , where the 2D bulk material would be non- 
topological. This resilience is due to the fact that the 
bulk p + ip and p — ip subbands are offset by different 
quantization energies. 

These results offer two principal advantages for realiz- 
ing Majorana fermions as end states in multi-band con- 
fined structures. First, the chemical potential in candi- 
date materials may naturally lie far away from \i = 0. 
For non-confined schemes, one would need to shift fi to 
lie in the Zeeman gap, possibly requiring chemical dop- 
ing or large electrostatic gate voltages. Here, one only 
needs to fine-tune the chemical potential on the order 

2 

of E s ]j ~ 2mL 2 ' Secondly, as discussed in Ref. 11, be- 
cause of the extremely small Rashba splittings available 
in semiconductor nanowires, restricting |/i| < V z in these 
materials would require operating at exceedingly small 
carrier density. However, the above results demonstrate 
that one may work far outside this |/i| < V z regime, open- 
ing the door to substantially higher densities. 

Despite this, small Rashba couplings still pose a serious 
problem for semiconductor materials. In order to remain 
in the topological regime, one needs V z > A. Given that 
typical values of ma 2 R in semiconductor materials are of 
the order of 0.2K 11 , and assuming that V z ~ A is on the 
order of a few Kelvin, one has V z /ma 2 R < 1/10. This 
means that only a small fraction of the superconducting 
gap would be converted into p-wave pairing, requiring 
one to operate at very low temperatures in order to avoid 
thermal excitations. 



V. FERROMAGNETIC DOMAINS 

We next consider a fully two-dimensional version of (7) 
in which V z is non-zero only inside a long finite strip of 
length L x and width L y . This could be accomplished by 
depositing a narrow strip of FM insulator on a fully 2D 
Rashba coupled surface. The idea here is that, under the 
right conditions, the ferromagnetic strip creates a strip of 
topological p+ip superconductor embedded in a topolog- 
ically trivially p±ip superconducting background. Such 
spatial interfaces between regions with different topolog- 
ical ordering generally give rise to localized zero-modes. 

Figure 6 shows a convenient geometry for realizing such 
a narrow domain wall. In this setup, the ferromagnetic 



insulator (FMI) layer has boundary conditions such that 
the magnetization points in the —y direction on the left- 
most side and in the +y direction on the rightmost side 
of the film. A narrow notch in the center of the film seeds 
the domain wall in the center of the sample. Since thin 
films of FMI are typically easy-plane materials (i.e. the 
magnetization prefers to lie in the x-y direction), it is 
energetically favorable for the domain to form by rotat- 
ing the magnetization out of the plane along the z-axis. 
Since only the z-direction magnetization opens up a Zee- 
man gap in the underlying Rashba coupled material, this 
domain wall creates a narrow topologically non-trivial 
region (shown as a shaded blue rectangle in panel b. of 
Figure 6). The resulting magnetization profile automati- 
cally varies smoothly with position, avoiding rough edges 
which would create pair breaking scattering. 

Fig. 7a. shows the bulk excitation gap as a function 
of chemical potential \i for an infinitely long FM strip 
of width L y = 20 embedded in a large non-FM coupled 
background with parameters t — 20, A = 10, olr = 10, 
and V z = 5. In this simulation V z has a sharp step profile, 
and drops abruptly from V z = 5 inside the strip to V z = 
outside. Red shaded regions indicate the presence of Ma- 
jorana end-states as determined using the Green's func- 
tion method (8). As with the case of hard- wall boundary 
conditions, we find an alternating sequence of topolog- 
ical phase transitions as \i is varied. These transitions 
again occur where the bulk gap for an infinitely long strip 
closes. However, unlike the hard- wall boundary case con- 
sidered above, k y is no longer quantized by confinement 
and the bulk gap closings are not determined by trans- 
verse momenta quantization. Rather, the spatially in- 
homogenous V z profile gives rise to an effective step po- 
tential which binds discretely spaced ID energy bands, 
whose bottoms lie within the bulk superconducting gap. 
Once again, Majorana end-states appear when an odd 
number of these discrete energy bands are occupied. 



a) 



x 



b) 



A 



FMI 



FIG. 6. (Color online) A smoothly varying ferromagnetic do- 
main wall can be created by creating a notch in a strip in 
an easy-plane ferromagnetic insulator (FMI), and seeding the 
opposite sides of the notch with opposite (in-plane) magne- 
tization. Panel b shows a top view of the proposed notch 
geometry, with arrows indicating the magnetization seed di- 
rections. Panel A shows a side view of the resulting magneti- 
zation profile (y- and z- components of the magnetization M) , 
which tips preferentially into the z-direction at the location 
of the notch, creating a thin region of non-vanishing z-axis 
polarization corresponding to the shaded blue region in b. 
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Fig. 7c.-e. show wavefunction profiles of the spin-up 
and down components of the Major ana end- states for M 
= 1,3, and 5 occupied subbands. The wave- functions are 
exponentially localized to the ends of the FM strip, with 
characteristic localization length set by the bulk excita- 
tion gap. The number of nodes in the end-state wave- 
function increases with J\f. Furthermore, since the V z 
potential step profile favors spin-down, the more deeply 
confined J\T = 1 and 3 states have spin-down, the shal- 
lowly confined J\f = 5 band has a more equal mixture of 
spin-up and down. 

In realistically fabricated structures where V z is in- 
duced by proximity to a FM insulator, the V z profile 
will likely not be a sharp step function. Fortunately, we 
find that the above picture is largely insensitive to the 
details Zeeman splitting spatial profile V z (y). For exam- 
ple, Fig. 7b. repeats Fig. 7a. with smoothly varying 
V z (y), obtained by applying a Gaussian filter with width 
of 3 lattice spacings to a sharp step-profile. The alternat- 
ing phase transition structure and presence of Major anas 
is very similar to the sharp step-profile case. The main 
difference being that the larger M sub-bands are more 
shallowly confined by the smoother V z potential. 

We have shown that the existence of Major ana end 




FIG. 7. (Color online) (Top) Bulk excitation spectrum for 
infinitely long ferromagnetic (FM) strip-domain of width L y 
with sharp edges (a) and Gaussian-smoothed edges (b). High- 
lighted red regions denote the presence of zero-energy end 
states as determined by the Green's function method (8) for 
strips of length L x = 100. Insets show V z (y) profiles. Tight- 
binding parameters used were t — 20, = 10, V z = 5, and 
L y = 20. (Bottom) Corresponding wave-function |^| 2 profiles 
for finite length (L x — 100) strips with sharp V z profile as in 
(a). Subfigures (c),(d), and (e) show |^| 2 for selected chemical 
potentials fi/t = —4.22, —4.12, and —3.87, which correspond 
toAf= 1,3 and 5 occupied sub-bands respectively. Spin-down 
(left) and spin-up (right) components shown separately. Faint 
dashed- lines mark the boundary of the FM strip. 



states in this setup is robust, persisting so long as the FM 
domains are narrow, elongated structures with gapped 
mid-sections. These results make it feasible to produce 
Major ana states in FM domains patterned on top of bulk 
2D Rashba split surfaces with induced superconductivity. 



VI. MAJORANA END-MODES WITHOUT 
TRANSVERSE SUB-BANDS 

A. Random Edge Geometries 

In the long rectangular strip geometries considered 
above, the system is neatly separable in the x and y 
directions. So far, the existence of discrete transverse 
sub-bands (in the y-direction) has played a central role 
in understanding the topological phase transitions in 
these structures. Naturally, one might therefore wonder 
whether the existence of transverse sub-bands is essential 
to the formation of Majorana end-states. Specifically, 
the presence of spatially varying and non-parallel edges 
mixes different transverse sub-bands, destroying the no- 
tion of the "number of occupied channels". Since, for 
rectangular samples, Majorana end-states exist only for 
an odd- number of transverse channels, it is possible that 
the mixing of even and odd number of channels may de- 
stroy the Majorana end-states in non-rectangular sam- 
ples. 

To address this question, we consider samples confined 
to a narrow region by electrostatic confinement potential 
Konf (x, y) with smooth, random meandering boundaries. 
Here it is important that the edge variation is relatively 
smooth, as jagged edge variations produce a scattering 
mean- free path t ~ W. Due to the condition W < £o> 
scattering from sharp edge variations tends to destroy 
the p-wave pairing gap 1 . 

To produce random edges with width variance aw 
and correlation length £, we start by choosing the y- 
location of the top and bottom edges yt,b(x) inde- 
pendently for each x, identically distributed normally 
with variance yla^l and mean L y /2 (where the over- 
bar indicates averaging with respect to edge configu- 
ration). We then apply an exponential smoothing fil- 
ter yt,b(x) — > Y2 X ' j e ~^ x ~ x ^ £ yt,b( x ')i which correlates 
yt,b( x ) an d yt,b( x ') on lengthscales, \x — x'\ < on the 
order of the edge-correlation length t. 

Since the sample width at any x must be an integer 
number of lattice spacings, {yt,b(x)} are rounded to the 
nearest integer, resulting in discrete steps rather than 
smooth edges. These steps introduce sharp, short range 
scattering potentials, and in order to separate out the 
effects of this discretization, from those of the smoothly 
wandering edges, we smooth the V con f along the lateral 
(y-) direction with a Gaussian filter of width 2 lattice 
spacings. 

Fig. 8a. and b. show the results of simulations with a 
smoothly random electrostatic confinement potential for 
samples with average width L y = 20, length L x = 100, 
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FIG. 8. (Color online) (a) and (b) show the two lowest energy 
in-gap excitations for an electrostatically confined strip with 
L x = 100, L y = 20, t = 10, A = 1, V z = 2, a R = 2. 
Red shading indicates the presence of isolated Major ana end- 
modes at zero-energy. The results in (a) are for straight edges 
(aw = 0), and those in (b) are for a random sample with 
aw = 4 and t — 15; (c) shows a colormap of the random edge 
geometry used to generate (b). Importantly the Majorana 
edge states survive, retaining a substantial excitation gap even 
for large edge variation (in this case ~ 40% of the average 
width L y ) and demonstrating that these states do not rely on 
the existence of transverse sub-bands. 




FIG. 9. (Color online) Spatial profile of the Majorana wave- 
function intensity, |^(x,^/)| 2 , for L-shaped junction, for t — 
10, A = 1, olr — 5, and V z — 2. The dimensions of each leg 
of the junction are 5 x 50 lattice sites. The left (a) shows the 
wave-function with three occupied sub-bands in each leg (/i = 
— 3t). In this case the Majorana states exist at the extreme 
ends of the the L. The right (b) shows the wave-function with 
three sub-bands occupied in the horizontal leg (/i = — 3t), and 
two sub-bands occupied in the vertical leg (/i = —3 At). In 
this case the second Majorana mode appears at the junction 
between topological and non-topological regions at the elbow 
of the L-junction. 



edge correlation length £ = 15, and with edge variance 
<Jw — and 4 respectively. We find that the Majo- 
rana end modes, and corresponding sequence of alternat- 
ing phase transitions survives even for substantial edge 
variations, that is, even when there are no well-defined 
transverse sub-bands. Despite the lack of transverse sub- 
bands, as one sweeps /i, discrete bulk levels inside the su- 
perconducting gap are still pulled down one-by-one across 
zero-energy, resulting again in an alternating sequence of 



topological phase transitions. However, the locations of 
these transitions occur at different values of /i compared 
to the rectangular case. In contrast to the rectangu- 
lar sample case, these discrete levels cannot be simply 
identified with transverse band-bottoms, but rather are 
bulk states with some more complicated structure. The 
excitation gap protecting Majorana end-states in these 
random edge geometries is reduced from the rectangular 
case. However, as seen by comparing Fig. 8a. and b., 
this excitation gap remains a substantial fraction of the 
straight edge rectangular case even for large variations in 
the edge geometry (in the case of Fig. 8 b. the fractional 
variation in width, (Sy^ + Sy^)/L y = 2aw / 'L y , is 40%). 

These simulations demonstrate that the existence Ma- 
jorana end-states is highly insensitive to the details of 
sample geometry, and in particular does not require the 
existence of transverse sub-bands. This robustness to 
edge- variations highlights the truly topological nature of 
these states. Also from a practical perspective, the abil- 
ity to tolerate substantial (smooth) edge variance eases 
the requirements for sample fabrication, making an ex- 
perimental realization more feasible. 



B. L-shaped Junction 

Majorana end-states can be moved, braided, and fused 
using appropriate networks of quantum wires 18 . A com- 
mon feature of such networks are right-angle junctions 
(either L- or T- shaped). Fig. 9 provides an explicit 
demonstration that the Majorana end-states with mulit- 
ply occupied sub-bands, even in geometries with sharp 
corners. Fig. 9 a) shows the wave-function for an odd- 
number of occupied sub-bands in both the horizontal and 
vertical legs of an L-junction, whereas Fig. 9 b) shows 
the wave-function for an odd-number of occupied sub- 
bands in the horizontal and leg and an even number in 
the vertical leg. In the former case, the entire L-junction 
is topological, and Majorana end-states reside at the ex- 
tremal ends, whereas in the latter case, the Majorans lie 
entirely within the horizontal leg, with the second Ma- 
jorana occuring at the junction between topological and 
non-topological regions. 



VII. CONCLUSION AND DISCUSSION 

In conclusion, by numerically diagonalizing the tight- 
binding model 7, we have shown that quasi-one- 
dimensional microstructures with Rashba spin-orbit cou- 
pling are are a robust medium in which to realize Majo- 
rana fermions. In spatially confined structures, such as 
semiconductor nanowires, confinement quantization en- 
ables one to operate at substantially higher carrier densi- 
ties where a bulk 2D material would be non-topological. 
Additionally, the presence of multiple energy scales al- 
lows one to trap Majorana fermions at the ends of quasi- 
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one-dimensional ferromagnetic domains, or in electro- 
statically confined strips. 

Importantly, Majorana end-states realized in this way 
are largely immune to both bulk disorder 1 and random 
sample geometry (so long as the edges are relatively 
smooth). Furthermore, for end-states, the mini-gap to 
localized excitation scales as A, rather than A 2 jep as for 
vortex core states. While it has been argued that mini- 
gaps are irrelevant for detection schemes based solely on 
electron number parity 14 (in which Major anas are kept 
far apart from each other), our results show that Majo- 



rana end-states allow for simpler fusion based measure- 
ment schemes. Finally, the ability to selectively tune seg- 
ments of these systems through a topological phase tran- 
sition, simply by electrostatic gating, provides a conve- 
nient, scalable avenue towards manipulating and braiding 
many Majorana fermions. 
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